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ABSTRACT 


Computer programs have been developed for calculation 
Of the molecular absorption and aerosol extinction for 
horizontal path propagation of radiation in the wavelength 
range above 0.55 um. The AFCRL line parameters compilation 
is used aS input data with local weather conditions for the 
Monterey Bay and model of typical aerosol size distribution 
From Shettle and Fenn. Wind speed and relative humidity 
are included in the form developed by Wells, Gal and Munn. 
Infinite resolution computations have been carried out for 
imvoo.c and 10.6 um and results tabulated for monthly 
and yearly average weather conditions on Monterey Bay. 
Bandwidth-averaged spectral transmittances have been calcu- 


~ aire On > com + resolution for the same 


ateca 2 0.01 cm 
SoncditLvoens. 

Comparisons with previously published data show agree- 
ment to 1.8% with maritime model aerosol Mie-scattering 
calculations by Selby using LOWTRAN III B [Ref. 12]. 

Analysis of sensitivity to changes of pressure and 
temperature for the 3.8 um range showed only small effects 
of these parameters. However a change in water vapor 


partial pressure results in a 0.75%/mbar change in trans- 


mittance in the 3.8 um region. 
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io LN eRODuUCT EON 


A. BACKGROUND 

aiemlight bringing Ws information from a distance, 
illuminating our surroundings and carrying energy, per- 
ZOerms these functions only by travelling in the atmos- 
PieGemeebaser Might an the form of a highly collimated 
beam with very narrow frequency bandwidth and high radiant 
power is equally dependent on the propagation properties 
Of the atmosphere. The characteristics of the light: 
motably its color (wavelength), intensity, polarization 
Miismopattat Glstribution are all modified by interaction 
with the constituents of the atmosphere - mostly the absorb- 
ing trace gases and suspended particles in the size range 
MemeOorene Order Of hundreds of micrometers. 

In modern military technology a wide range of systems 
are in use or under development for surveillance, imaging, 
eameee Getection, designation, identification and tracking, 
PeEwecetmmunicatlons, for precision guidance of munitions 
and possibly even as a directed energy weapon. These oper- 
ate over the entire "optical" range of wavelengths, from 
the visible or near ultra-violet through the "far infra red" 
from about 0.3 to 20 micrometers. Operation of these sys- 
tems is limited by the properties of the atmosphere as a 
transmission medium, particularly the degradation of 


resolution by scattering from particulates and turbulence, 
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and loss of radiant intensity through absorption and scatter- 
ice DUGEHto strong absorption by atmospheric molecules, 
Systems are constrained to operate in limited wavelength 
ranges or "windows", relatively free of absorption lines. 
Of these windows, the radiating temperatures of many 
sources of military interest and the desirability of day/ 
night operating capability have concentrated attention on 
the 3-5 um and 8-14 um infra-red windows. 
[Mmraretcular the navy has a special interest in the 
DF iaser aS a candidate moderate or high power laser at a 
Wavelength in the 3-5 um range, unaffected by the high 
humidity levels found in the naval environment. Of particu- 
lar interest is the 3.8007 P. (8) EGanoile Om lines wiilcl 
Shows up to 96% transmission over a kilometer path at sea 
level. 
The need to predict the operational performance of 
these systems leads to a requirement to predict the spec- 
tral transmittance of the atmosphere in the infra-red, as 
a function of the meteorological parameters which determine 
the absorption and scattering levels — the pressure, tem- 
perature, air composition, vertical gradients, number den- 
sity composition and size distribution of suspended particles. 
Large scale meaSurement programs have been undertaken 
in many countries with the objective of building a large 
miprary Of Meteorological and transmission data. Consider- 


able effort (notably at AFGL) has been devoted to the 


cea 





development of computer programs for transmittance calcu- 
lation based on weather data, and the establishment of a 
number of "standard" atmospheric models. 

The current modelling effort was undertaken in support 
of the experimental program of transmittance measurement 
on ranges over Monterey Bay being carried out by the 
Environmental Physics Group at the Naval Postgraduate 
school. This work involves measurement in all the suitable 
atmospheric windows from the visible to 14 um, with both 
laser and broadband sources. The calculational program 
has therefore been directed to computation of molecular 
absorption and aerosol extinction coefficients for sea 
level propagation under a range of weather conditions sta- 
meoetecaliv typical Of Monterey» This has been done for a 
variety of wavelengths appropriate to the measurements, 
and at high resolution for use with laser sources, and 
averaged over filter bandwidths for broadband sources. The 
Air Force Cambridge Research Laboratory Atmospheric Absorp- 
tion Line Parameters Compilation [Ref. 7] in magnetic tape 
form is used as input data for the absorption calculations. 
A Mie scattering program based on the Van De Hulst approxi- 
mation with the model particle distribution of Shettle and 
Fenn [Ref. 8] has been used for the aerosol scattering. 
Wind speed and relative humidity effects have been included 
uSing the drop size growth model of Wells, Gal and Munn 
[Ref. 18]. A sensitivity analysis has been carried out with 


respect to pressure and temperature. 


eZ 





Basically, this thesis describes a program to 

a) Develop a computer code to calculate atmospheric 
transmittance for specific wavelengths in the near and 
mad IR. 

6b) Develop a code to give transmittance and molecular 
absorptance with bandwidth appropriate to grey-body sources 
in the near and mid IR. 

emeopely these £o prediction of transmittance over 
marine optical paths, over Monterey Bay, based on the 
available meteorological data base for the area, for con- 
ditions appropriate to the optical propagation experimental 
MeasuLement program. 

Curves have been prepared giving predicted molecular 
and aerosol extinctions for statistical monthly and yearly 
average meteorological conditions on Monterey Bay. 

Further sections of this thesis describe the theoreti- 
cal background and the details of the algorithms, and provide 
Eebular and graphical results of the calculations. Computer 


program listings are given in the appendices. 


B. ATMOSPHERIC PHENOMENA 
Optical transmission of energy is usually described 


by the transmittance, 





tT = exp(-uL) = (1-1) 


ILS 





where u is called the extinction SseGereient, and wil yepre— 
Sents the length of the path which the light has traveled. 
This expression follows from Beer's Law, which for linear 


propagation of monochromatic radiation at frequency v has 


eae. form 


= GZ i) G7) a) 


Here I(z,u) represents the intensity of monochromatic 
meelatton at a distance z. 

Beer’s Law in differential form (1-2) may also be 
written as 


niu =o tos eee G3) 


BPeaemaegeven Path length L-. The transmittance Tt over a 


path of length L is then given from the equation (1-3) 


Wes Gi == Se Me 
TU), = TE 2 


which is the same as equation (1-1) for monochromatic 
maclat i On . 

The extinction coefficient »# is generally the result 
Of two processes: absorption and scattering. Both mole- 
cules and particles suspended in the atmosphere cause 
absorption and scattering. It may be useful to define the 


SsEimetlon cCoerrLicient u aS a sum of four terms 
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u = k + g + K + 2 (1-4) 


where Ka! Ba! Kae 8. are molecular absorption, molecular 
scattering, aerosol absorption, aerosol scattering 
Soctiaecilents. 

An aerosol extinction coefficient Ue 1s sometimes used 
for the summation of aerosol absorption and scattering. 
eemebarcty a molecular extinction coefficient us is used 
for the summation of molecular absorption and scattering. 


From equations (1-1) and (1-4) the total transmittance 


through a path in the atmosphere may be written as 


T({total) = at(molecular absorption) * 
Mine mee wlan Scaecerl ig) =< 
T(aerosol absorption) ~x 


t(aerosol scattering) 


where t(molecular absorption) may be written in terms of 


the atmospheric components as 


t(molecular absorption) = t(uniformly mixed gases line abs.) 
TLOZOne eine abs.) < 
t(water vapor line abs.) * 


Ti water vaoon Con-Elnuum abs -) 


Pineesoqemmceneltmuum abSay) 


ILS 


x 
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Although the molecular Seaereming ~ederosol abscorae lon 
and aerosol scattering coefficients are smooth functions 
of frequency of monochromatic radiation, the molecular 
absorption coefficient is strongly dependent on frequency. 

Molecules which cause absorption and scattering of a 
light beam in the atmosphere may be divided into two types: 
uniformly and non-uniformly mixed gases. Carbon dioxide, 
nitrous dioxide, carbon monoxide, methane, oxygen and 
nitrogen are assumed as uniformly mixed gases for which 
their Mixing ratios by volume with each other are constant, 
independent of pressure and temperature. These mixing ratios 
are tabulated in table 4. Water vapor and ozone are con- 
Sidered as non-uniformly mixed gases, whose relative con- 
centrations must be specified. 

Six different models of the atmosphere are frequently 
used in calculations; these are the tropical, midlatitude 
Summer, midlatitude winter, subarctic summer, subarctic 
Winter and U.S. Standard atmosphere 1962 model described by 
McClatchey and D'agati [Ref. 4]. Pressure, temperature, 
density, water vapor and ozone concentrations in the air 
for each model at sea level are tabulated in table 5. 

The number densities of aerosol particles which also 
cause absorption and scattering of a light beam are more 
complicated to determine. It is necessary to know the size 
distribution of aerosol particles in addition to the number 


density, in order to compute the aerosol scattering and 


kG 





absorption coefficients. The size distribution function 


is denoted by n(r) and given as 


Be) Gale) 


where N(r) is called the number density which is the number 
of aerosol particles with radius less than r perunit volume. 

Several different analytical models have been developed 
meee Size distribution. A “continental” aerosol model, 
Whteh includes both rural and urban aerosol conditions, 
and a "Maritime" aerosol model are used for the lower 
atmosphere. 

The rural aerosol model used by Shettle and Fenn is 
assumed to be composed of a mixture of 70% water soluble 
Substance, consisting of ammonium and calcium sulfate and 
organic compounds, and 30% dustlike aerosols [Ref. 8]. 

In urban areas the rural aerosol background is modified 
by the addition of aerosols from combustion products and 
industrial sources. The proportions of soot-like aerosols 
and the rural aerosol type aerosol mixture are assumed to 
be 35% and 65% respectively [Ref. 8]. 

The maritime aerosol model represents the aerosol com- 
position and size and number distributions over the oceans. 
Peioworomttlcantly different from the continental aerosol 
types. The maritime model is due to salt particles caused 


by the evaporation of the seaspray droplets also added to 


ey 





the rural type aerosols, from which most of the very large 
particles have been eliminated. 
The general form of the size distribution for these 


Mmeaels 1S given by Shettle and Fenn [Ref. 8] as 


; (Eege iT log BG 
1 ee = 
Beat ira 10) ro, (27) 20; 
(1-6a) 
where the parameters in the equation above are tabulated 
in table 6. 

The size distribution given by equation (1l-6a) is 
relative humidity and wind speed independent. Since the 
number of oceanic particles with greater radius in the size 
distribution increases with increaSing relative humidity and 
Wind speed, a modified maritime model has been proposed 
Dmom@eecadleulation of aerosol extinction coefficients taking 
SamemOor the growth Of particle radius {Ref. 18]. This 


gives 
me K(0.47 (5) ~* + 2.3Mpa(ar)exp[-8.5(er) 1} (1-6b) 


emere CC 1S a normalization constant, and Mp ls the mixing 
ratio of seaspray aerosols with the rural aerosols. Or 
Poememeteathy tne Mixing ratio of rural-Oceanic aerosols 
PimeteweauacLlon (1-6b) is 1:M,- 


Particle growth is included through the growth factor ¢ 


which adds the effect of relative humidity into the equation 


lle: 





@eago). lt 1S given as 


Pe ee Oo, GO = 





100” 


where R.H. 1S the percent relative humidity. The effect 


of F 1S to increase the drop radius to 


where re 1S the radius for zero percent relative humidity. 
F iS normalized to 80% relative humidity, given F'. There- 


menbenwtactcor F*' in the equation (1-6b) is 


ECR. Aw 
F (80) 


Y in equation (1-6b) describes the dependence on the wind 
Speed. It has the form 


Me TD, Bh MO OURS aa oe 


where V is the wind speed in meters per second. ain 
equation (1-6b) also represents the wind speed dependence 
and it iS given as 


ae oO. te 7 50 yr? 16 ; Wien 7 imn/ SEC. 


Bo to when V > 7 m/sec. 


ow 
lI 
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So the inputs for equation (l1-6b) are wind Speed V and 
relative humidity R.H. The normalization constant C may 
be determined from the relation between visibility and 


aerosol extinction coefficient at 0.55 micrometers. 


C. HF/DF LASER AT THE NAVAL POSTGRADUATE SCHOOL 

Although this thesis will deal with propagation at 
Pea level at any frequency of radiation, primary interest 
will be the frequencies of the HF/DF laser. These fre- 
quencies are tabulated in tables 7 and 8. 

To be used for the propagation measurement program at 
the Naval Postgraduate School, a small scale hydrogen 
Fluoride (HF) or deuterium fluoride (DF) laser was located 
Sieene FOOL Of Spanagel Hall. The system consists of a 
vacuum SyYStem, a power supply and the laser head. The 
laser is centered around a laser head which contains dis- 
charge tube, mixing chamber or cavity and resonant cavity. 
The system uses four gases: a) Sulfur hexafluoride, b) oxy- 
gen, c) helium, d) hydrogen or deuterium. A mixture of 
sulfur hexafluoride, oxygen and helium is injected into the 
discharge tube at the end opposite the mixing chamber. 

Then the sulfur hexafluoride is dissociated by electron 
impact in the discharge tube by making use of an electric 
discharge. Oxygen helps the sulfur hexafluoride to dis- 


sociate into free fluorine radicals by the reaction 
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The hydrogen or deuterium is injected perpendicular to the 
flow in the mixing chamber. Then chemical reaction takes 
place in this mixing region forming the vibrationally 
excited HF or DF molecules. 

The optical resonator of the system consists of a 
pair of Brewster windows and mirrors. The Brewster windows 
which linearly polarize the output are chosen from two inter- 
changeable pairs, one of BaF. (025 — 15m ecransmieci ng) 
piceene OLher sapphire (0.14-6.5 um transmitting). 

The system has two options for operation: a) Multi- 
ime @oetation, b) Single line Operation. The output 
mieror should be replaced by the diffraction grating for 
Single line operation, while two mirrors (one partially 
transmitting) are used for multiline operation. 

Optimized power output from lasing of hydrogen 
fluoride was reported [Ref. 9] as 6.5 Watts with the mass 


flow rates 


—" 
try 
i 


0.600 grams/sec. 


6 
0, = 0.182 grams/sec. 
H. Oe Oona teams) Sec. 
He SOR Oa Ge git ams, see. 


Optimized multiline power output of deuterium fluoride 
lasing was also reported [Ref. 9] as 2.5 Watts with the 


mass flow rates 
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SF = 0.680 grams/sec. 


6 
0. = 0.270 grams/sec. 
D. = 0.043 grams/sec. 
He = 0.063 grams/sec. 


In the case of single line operation, the wavelengths 
and corresponding relative intensities of fourteen individual 
DF and fifteen individual HF lasing transitions are shown 
mmetgure 3 and figure 4. 

The block diagram of the main elements of the system 


1s also shown as figure 5. 
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II. OPTICAL PROPAGATION THROUGH THE ATMOSPHERE 


A. MOLECULAR SCATTERING COEFFICIENT 

The simple Lorentz-type molecular model presented here 
is shown in figure 1. It is assumed nonionized, nonpolar, 
1sotropic, linear and lightly damped. 

The model contains molecular mass residing in the 
atomic nuclei and carrying a net positive charge. Orbital 
electrons shared by the constituent atoms are represented 
by a concentric shell, providing a balancing negative 
Charge. Force directed toward the central mass is char- 
acterized by a spring constant k. Let m be the electron 
mass. If the system is disturbed by an electric field E 
Gn the z direction, then the equation of motion may be 


written for a displaced electron as 





d“z 
= cE, - kz = m 5 (2-1) 
ade 
where 
E =E Sin wt. 
Zz Sz. 
Beeeton (2-1) can also be written as 
2 =eE 
Se wiz = << Sin wt (2-2) 
elie 
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e/a 
where i (k/m) / 1s defined as the resonance frequency, 
corresponding to a molecular transition. 


iiewsteady state solution to equation (2-2) is 


as 


If we denote the maximum value of the dipole moment 


aS Por emen 1 1s 


The oscillating dipole moment has maximum value Po 
given by equation (2-3) and produces secondary waves. 
Reradiation from the oscillating dipole then produces 
Secreto and Magnetic fields at every point in the coor- 
dinate system shown in figure 2. In the short dipole 
approximation (i.e., radius of scatterer is much less than 
wavelength of incident wave) and far field (i.e., distance 
to scatterer is much larger than radius of scatterer) the 


electric and magnetic field vectors can be written as 
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A O Sin 6 .. x 
le = an = Sin (kr-wt) 8 
> ao 

= Oo See : a 
B lg errr Sin(kr-wt) 6 


where 


ers Liem elOCLEy OLewave, 
HL is the magnetic permeability, 


8, ¢and r are the variables in spherical coordinates. 


‘ ae s 
The poynting vector denoted by S describes scattered power 


per unit area, and is given here as 


2 
4 EXB ye), 2 cea 8, ARN? 


S = ————— = ) ame Sin 


a ae Chui ae (2-4) 


Scattered power per unlit solid angle can be found from 
pPem@ecton (2-4), first taking the tlme average, then multi- 


ping by -. This gives the scattered intensity as 


2 
WL BG ee 2 


u WwW e 
a ( )"E 
ev 40 4mm (we - a) 





(2-5) 
the power of the incident wave contained in unit area can 


be given as 


=a Vanes ee (2-6) 


a ; 
<s incident 2 Oz 


where € is the permittivity of the medium. 
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The angular cross section of a molecule, denoted by 
gf 8)s may be defined as that cross section of incident 
wave acted on by the molecule, having an area such that 


Poevwe= tTlowing across it is equal to the power scattered 


Myeeece molecule per unit solid angle at an angle §. It 
1s 
no) = Does?) fee esac ton 2) ae 
m Sole nana equation, (2-6) 
Z 2 Zz p 
Uae 2 Z 2 
a8) = (Ne = EES (2-7) 
: (ee ) 2 


It is necessary to define the polarization P in order to 
find a suitable form for the second term in the right hand 
Side of equation (2-7). P is defined as total dipole 


moment per unit volume and it is given as 


for a linear and isotropic medium, where 


N is number of molecules per unit volume, 
ye is electric susceptibility, 
K 1s dielectric constant, 
EOFF iowa efrective Glecernte Fiera which is 
Eon * Jer 
7 O 
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Now we can write P in terms of K and Eo , as 
Z 


3(K-l)e. 
a K+2 Bas i 
fcom woiich 
ie _ 3(K-l)e, 
E N (K+2) 
OZ 


The index of refraction is given here as 


ee ee 4 ee 
=) KK VK 
Since Kn = 1 for a nonmagnetic medium. 


Now equation (2-7) can be written as 


Sea defined by the equation below, is called the total 


Scattering cross section. 


an 
oe tera) cle 
0 
OF 
3 2 
go = (stg)? (2-8) 
= Nee lee) 
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where 


was substituted. 

The total scattering cross section given by equation 
Wee) iS defined similarly to the angular scattering cross 
section, as the power scattered in all directions. 

For air the refractive index is very close to one. 
Therefore meee is nearly equal to three, and equation (2-8) 


may be written as 
—>=—, (n -l) (2-9) 


Mees £Orm of the total scattering cross section (2-9) 

has been derived for a single molecule disturbed by incoming 
Pecenomagnetic radiation. It has the units of area. For 
many molecules in a unit volume, the total scattering 
Seoerficient =n may be defined, assuming each molecule has 


the same total scattering cross section, 
ee 
B =. Neeee= ———- (n--1) (2-10) 


Hiemmolecular scattering coefficient (2-10) has the units 
eel /length. 

Since atmospheric gases are not quite isotropic, a 
@eeeolari zation factor Ps is introauced (Ret. 1] and inserted 


mieomeguation (2-10). The molecular scattering coefficient 


28 





may then be written as 


m 4 6 7P Pp ik 
3NA O 


taking into account the pressure and temperature dependence. 
Here aa and or refer to standard conditions, 1013 mb, 
fo Wk. 

Assuming air to be an ideal gas and using the ideal 
gas law, air is found to have a mole density of 44.631 moles 
per cubic meter at a pressure of 1913 millibars and tempera- 
ture of 273 degrees Kelvin. This is equivalent to 
2.688 x lige molecules per cubic meter. 

For a depolarization factor, a value of 0.035 is 
used, following Penndorf [Ref. 1]. Although the index of 
refraction of air depends on the wavelength, this variation 
is small. Therefore eye 1s taken as 3.0884 x ti 
which corresponds to a wavelength 0.55 micrometers [Ref. 2]. 


Using these parameter values given above, the atmospheric 


molecular scattering coefficient may be written as 


a 1OROmaS a. P 
om 4 1013 
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where 


MaSceum:ts Of “cm. 


>? 


P has units of mbars. 


T has units of degrees Kelvin. 
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McClatchey and d'Agati [Ref. 4] have used the slightly 


@eererent form of the molecular SCateerIng cocer1e1ent, 


(2a) 


which was selected as giving a "best fit" to the experi- 


mental data of Penndorf [Ref. 6]. 


Bee MOLECULAR ABSORPTION COEFFICIENT 
The model to derive the molecular absorption coeffi- 

cient is the same as the one we developed in the previous 
section, but now it is heavily damped. When the model is 
disturbed by an electric field, the electron moves with a 
@ampeinig coefficient ance restoring EOrce Strengen rk. The 
equation of motion may be written in differential form such 
that 

a*r ar 


> 
m—oommt+my = + kr = -eE = -eE(r)e 
dt ae 


NO 


Miewsteady state solution to this equation is 


= - 
2 so | 
woot we 7 Lye 


_> 
E 


The polarization due to N molecules per unit volume is 


Nh 
¥ 


p = —-eNr = — = BE Zo) 
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where we again replaced the applied electric field with 
an effective electric field to account for POldrizaotl tty 


of the medium. 








Meeer SUbStituting E = E + P/3e , the polarization 
eff O 
becomes 
Se oO ae 
Pp = 5 E 
Wy, - wo = lyw 
where 
Moose e“N 
1. o 3e 5m 


Pim~eroe POlarization is also given as 





= > ~ 2 > 
WE on ae 8 ere = Ag ee ae 
(2-13) 
Semmiming equation (2-12) and (2-13) one can get 
Spe m 
a OP lea me aoe O 72 
1 aol Sb ear ere 
e) ae os 2a = ule gee 
(2-14) 
Bomenee Indaex Of refraction. 
Eemene general form for equation (2-14) is written 
below with j referring to different types of molecules 
Weeem Oscillator strength fs. 
2 N. £- 
J 2 
; ote aa ya (Desp 
€_m Z ; 
O W = 70) ey -U 
; iL 
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Near a resonance frequency, equation (2-15) can be 


approximated by 


using the Binomial Expansion, Since n is very close to 
Gaity . 
A plane wave going in the z direction with wave vector 


K may be written as 
—- ~~ 
i eee 


(Dae) 


The relation between wave vector and refractive index 


iS 
ot) 
S we 2-18 
K pod ( ) 


K is complex, since n is complex, so defining K and n 


A 
Il 


8 + 10 (2215) 


n + ik (2-20) 


Zz 
ll 


where 8, a, n, k are the so-called propagation constant, 
absorption coefficient, index of refraction, and extinction 


coefficient respectively. 


By 





Using equations (2-16,18,19,20), it can be shown that 


7 Nyf, wr ¥ 
a = (2=7 4) 
Ze eer (W547W ) sei 7 


Substituting K from equation (2-19) into equation (2-17) 


one can get 


bly 
Hl 


> ah . = 
Eo 2 AZ ai (sz tite.) 


The electric field decreases by a factor of Barer ag 


and the intensity falls off by a factor of ee 


along the 
path length z. Since a corresponds to absorption, many 
workers use the so-called Lorentzian line shape given in 


equation (2-21) in order to find an expression for molecular 


mesOmpelon coefficient. 


Substituting the expressions for frequency w = 2tTv and 
45 = 20V 5 into equation (2-21) a becomes 
A 
a = 
Aq 
=—- - + 
ie Oe] 1 
or 
2 
— A Omens (2-22) 


33 





where A is a constant, Av is y/4m and defined as half 
width at half maximum. 

Ifa is normalized in the range from minus infinity to 
plus infinity, the expression for the Lorentzian line 


shape becomes 





fee Deng Ene Spectral distribution of absorption. 
mercy tie molecular absorption coefficient can be 


written as 


ee 
TG) eo] 
Peovided S = JK (v) dv for single absorbing line, where 


Sues line eee” aay per absorbing molecule in 
eae UnLeS Of crit /(molecules/cm?) , 


Wet sSeoONNe!: width at halt maximum in the units 


Ot “CmEr- ; 
ee : ail 
yY 1S incoming wave frequency in the units of cm ", 
ve is resonance frequency ben waaseocping line In 


the units of cm 


The radiant frequency v 1S unique: this gives absorption 
as a spectral quantity. Radiation from laser sources has 
very narrow line width, which may be assumed infinitely 

narrow for calculation purposes. The resonance frequency 


for each absorbing line is independent of pressure and 
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temperature. The line width is assumed to be a pressure 
broadened Lorentzian, with its pressure and temperature 
dependence given by the expression [Ref. 7] 

al 


a = a(Po,T.) ge (ae)? (eos 
O 


where 


rd 
i 


1 atm. and 


296 degrees Kelvin. 


| 
ll 


The total intensity of any line is pressure independent 
and its temperature dependence is given by the expression 


[Ref. 7] 





fb 
S(T) = S(T.) nS sss x exp[1.439E" ( 


where Qy and Oe are vibrational and rotational partition 
functions respectively, E" is the lower state energy of 
mien ecansi clon corresponding to that line, in units of ae 


Rotational partition functions are also temperature 


depenaent, of the form 


The corresponding values of Qn and Geese tabulated for each 


molecule in the air [Ref. 7] and presented in table 2. 
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Finally the total molecular absorption coefficient due 
to water vapor, ozone, carbon dioxide, methane, carbon 
monoxide, oxygen, and nitrous oxide in the air is given 


by 


5,54, .m. 


Summation over j} corresponds to different types of mole- 
Sutes, and Summation over 1 corresponds to different lines 
for a given molecule. 

Values of Ug a ay or ee S54 (T.) have been compiled 
for 100,000 lines between 1 um and the far infrared by 
McClatchey and others at Air Force Cambridge Research 
Laboratories and are available on magnetic tape [Ref. 7]. 
Also a new version of the compilation including over 
139,000 lines between 0.2 um and 30 um was recently reported 
pyeeae FOrce Geophysics Laboratory [{Ref. 3]. 

Saenlation for i at any pressure and temperature may 
be made by using equations (2-23,24,25) with the data from 
this tape. 

m. in equation (2-25) represents the number of j type 
molecules per unit volume. The resultant K is in the windts 
of l/length and generally 1/kilometer is used. 

It should be noted that the Lorentzian line shape must 


be modified for pressures less than 100 millibars; at pressures 


between10 and 100 millibars doppler broadening becomes 
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dominant [Ref. 7], while at pressures less than 10 milli- 
bars, the Voigt profile should be used. The Lorentzian 
line shape is good enough to be used at low Gee UG es , otic 
line wings should be truncated at 20 om from the line 
centers, Since uncertainties begin to arise beyond that 
range. 

A computer program calculating the molecular absorption 
coefficient is prepared and presented as Appendix A. This 
program also includes the water vapor continuum absorption 
Gaiculation which will be explained in the next section. 
Mee aig@rithm for the program is also discussed in Section 


Jie 


Ce CONTINUUM ABSORPTION COEFFICIENT 
Continuous regions of absorption in the atmospheric 


t (C4 Eibaa)) a veligtel) aesereygg) 7 18)(0 ai EO 


windows occur near 2500 cm. 
1250 cm > (8-14 um). 
Iee@rtker Spectral regions, the contribution from nearby 
eesorotion lines is much greater than that from the continuum 
absorption, so that for practical purposes the continuum 
effect can be neglected, although it may be greater than 
in the windows. 
Extreme wings of strong collision broadened absorption 
.lines centered more than 20 em > away or pressure induced 


absorption resulting from transitions that are forbidden for 


unperturbed molecules may cause continuous absorption. In 
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mae oO-l4 um region, existence of the water dimer (H,O:H,0) 
Mey also cause continuous absorption. 
1. Water Vapor Continuum 
The absorption coefficient due to the water vapor 


continuum can be described empirically by the expression 


[Ref. 4] 


vas 
= 
i 


[Co (v,T)P. + Cy (v,T) PW (2-26) 


where Co 1s a self broadening coefficient due to self 
broadening from collisions of water molecules with other 
water molecules, Cy ese COre ton sgas beOadenIng jocrtieient 
due to nitrogen broadening from collisions of water mole- 
cules with air molecules which are mostly nitrogen. 

Co and or are both temperature and frequency dependent. 

epemo-14 wm region [Ref. 4] 
Gms lcion spy Ene expression 


N 


Cy (v,T) = 0.002 x Co (v,296). 


oe is given by the expression 


ee 


796) | 


Dee = Casa 2.o Oe exp [1800 (= ~ 


where Colv ,296) is given as 
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C.(v,296) = 4.18 + 5578 x exp(-7.87 x 107? v) 
where Co 1S in the units of (precipitable-cm) ~) ene or 
@Pr. aie ~ stm ~ which is (grams/cm*) ~1 een. 

Eee. 5-4.2eim region [Refyod] 


C, and Cy are given by the expressions below; 


Cy (v7 T) = 0.12C.,(v,T) 
Cet)! =. CH ,296) 9x exp [1350 (= 2 reas 
Saad Sie Sle 256 
where Co mas the units o£ (pr. em) 7+ atm >. The values of 


Co (v,296) at different frequencies are reproduced [Ref. 4] 
and tabulated in table 3. W in the equation (2-26) has 
the units of (pr. cm) per kilometer which is the same as 
grams per square centimer per kilometer. Then the water 
Vapor continuum absorption coefficient K can be 


-1 GOner 
Obtained in (kilometers) 


2. Nitrogen Continuum 
Although the nitrogen continuum is not included in 
this thesis, the magnitude of this effect is of the order 
Dyin AO ee Oval km > at a wavelength of 3.8 um. The maximum 
Setect walen 1S °O0f the order of O.1 em > occurs at a wave- 


length of 4.3 um. This effect may be neglected outside 


mee se7—2¢.7 um region [Ref. 4]. 
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Dee SEROSOL EXTINCTION COEFFICIENT 

Aerosol extinction includes absorption and scattering 
Byeeeecrosol particles. Although scattering by aerosol 
particles whose radius is smaller than about 0.03 times the 
wavelength of the light may be calculated by Rayleigh 
theory, for the whole range of particle radius Mie theory 
must be applied. The bases for Mie scattering are the 
interactions between electromagnetic waves and the elec- 
Gime charges that constitute matter. In the molecular 
scattering case which was described before, only a single 
dipole is involved, but in the case of aerosol scattering, 
a particle consists of many closely packed, complex mole- 
cules. Thus it may be considered an array of multipoles, 
Since an aerosol particle is a dispersed system of small 
particles suspended in the air, such as combustion pro- 
ducts, dust grains, bits of sea salt surrounded by water 
mebeemres, volcanic ash, etc. The multipoles give rise to 
secondary electric and magnetic partial waves. In the far 
field two kinds of waves combine to produce the scattered 
wave. Interference between two waves is mainly determined 
bpeeene wavelength Of incident light, and the size and rela- 
mEyecwGerLbactive index of the particle. A useful parameter 
used in Mie theory is the size parameter a, a dimension- 


less quantity defined by the relation 


SSeS IES (22975) 
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Mm@eme es 15 the radius of the scatterer or absorber, and i 
is the wavelength of the incident radiation. a is the 
parameter on which aerosol scattering and absorption depend 
Mmese sclOngly, rather than on radius of particle alone or 
6m the wavelength of incident radiation alone. The refrac- 
tive indices used in calculations are the relative refrac- 
tive indices of the aerosol particles with respect to the 
refractive index of the medium surrounding the particles. 
Since air has refractive index very close to unity, the 
relative refractive index is assumed equal to the absolute 
refractive index of the aerosol particle, and is generally 
represented by a complex number, where the imaginary part 
mines cate Of possible absorption by the particle. The 
refractive indices of different kinds of particles used in 
Calculations are reproduced and presented in table 10 
(Ref. 12], taking into account the wavelength dependence. 
Tercrleticght scattered of scattered and absorbed by a 
pig le particle is given by its total scattering or 
extinction cross section which is related to the geometric 
cross section by a parameter generally called the effi- 
Bileticy factor. According to the Mie theory, the scattering 


and extinction efficiency factors are given by the equations 





"$0 2 2, 
eS ) (2n+1) []a,| *+|b_ | (2-28) 
et & Bs eI 
ex 2 . , 
nea = y= =a ) (2n+1) [Re (a,+b.) | (2-29) 
i es 


n=l 





iene efficiency factor is defined as the ratio of 
apparent cross section to geometric cross section of the 
mactcicle. i eee in the above are the apparent scattering 
Emig extinction total cross sections respectively. Re 
represents the real part of the argument. The coefficients 
ane bo ane the complex functions which represent the 
amplitude of sae electric partial wave and ae magnetic 
partial wave respectively. These coefficients may be 

Beume from Ricatti-Bessel functions for a given complex 
ioiaadeceive index and size parameter a. 

The scattering efficiency factor approaches the value 
of two in the short wavelength approximation, which implies 
that the scattering cross section is almost twice as large 
as the geometric cross section of the particle. 

Since the infinite series represented by equations 
Meese and (2-29) are slowly converging, calculations are 
generally computer time consuming. For particles with real 
part of refractive index between 1 and 1.5 and imaginary 
part of refractive index between 0 and 0.25, the Van de 
Hulst approximation may be used with some correction fac- 
tors. Within the approximation range, it was reported by 
Begemenaygitan (Ret. 13] that the efficiency factor could 
be predicted within the precision range of + 0.05 oS: 

The method of calculation using the Van de Hulst 
approximation is presented in Appendix C. 

The efficiency factor of a single particle then may 


be extended for uSe in calculation of the aerosol extinction 
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coefficient due to many particles for a given size 


Bese rtbution such that 


nD 
8. = ze (oe Q.. nla) da B20) 
ail 
“D 
u, 7) f a* Q. nla) da (Q=an) 
al 


where K 1s the propagation constant and n(a) is the size 
distribution defined by equations (1-6) and (2-27). 

Once the aerosol scattering and extinction coefficients 
are determined, then the aerosol absorption coefficient may 
wemround from Se ae = 3a if needed. 

Size distributions used here for calculation of aerosol 
extinction coefficient are a) equation (1-6a) which is 
independent of wind velocity and relative humidity having 
75% sea spray and 25% rural aerosol mixture which is com- 
posed of 30% dustlike and 70% water soluble substance. 

b) equation (1-6b) which adds the effects of wind velocity 
Ee@wrelative humidity. 

If the total number of particles 1s not known, then the 
Size distribution n(a) or n(r) 1S generally normalized to 
one particle per cubic centimeter. Whatever value the 
size distribution is normalized to we can put a normaliza- 
MEOWmcenmstant into equations (2-30) and (2-31). Thus we 


may write 
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a S me Joa 2 igi eR) weteke! 


where C is the normalization constant and may be found 


meom the relation [Ref. 19] 


59 
B (0.55 um ) 


Pere R is the visibility in kilometers providing 8 (0.55) 
in inverse kilometers. 

Therefore knowing the range R for visibility, the aero- 
sol scattering coefficient can be found from the equation 
pgeve-. Imac acrosol scattering coefficient at 0.55 micro- 
meters calculated by equation (2-30) then allows determination 
Of the value of the normalization constant C. Once C is 
determined at 0.55 micrometers, it is the same for the 
whole range of wavelengths. 

Since aerosol absorption by particles is much less than 


Scattering, the visibility range may also be written as 


3.9 
We (Css ein) 


The curves of extinction coefficient versus wavelength 
for varied wind speed and relative humidity shown as figures 
imagine ace nmormalized so that the total aerosol extinction 

i 


coefficient at 0.55 um is 0.17 km. This corresponds to 


emvisteility of 23 kilometers. 
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ie aeCOR DTA OR CALCULATION -OF 
MOLECULAR ABSORPTION COEFFICIENT 

The computer program presented as Appendix A is 
intended to calculate the molecular absorption coefficient 
including the water vapor continuum effect. As described 
in Section II.B, the program reads the line parameters 
as? ee Seay from a magnetic tape. However, other 
inputs of the program are the amounts of the air molecules 
which have to be specified in the units of molecules per 
square centimeter per kilometer path, and they should be 
Supplied by the user. 

The amounts of different types of air molecules excluding 
the ozone and water vapor are calculated from the ideal gas 


law, 


where P is air pressure in Newtons per square meter, I 

is air temperature in degrees Kelvin, R is the gas con- 
Stant for air in Joules per kilogram per degree Kelvin 

and has the value of 287.06 Joules per kilogram per degree 
Kelvin. For given values of pressure and temperature, the 
density of air can be found in the units of kilograms per 
cubic meter. Then the density of air is used to find the 
amounts of the uniformly mixed gases (O., No, CO,, Cor 

N 


eeeen) such that 


2 + 
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where X, 1s the concentration in PPM by volume taken from 
table 4, and Os 1s the density or the mass per unit volume 
of the corresponding constituent. After finding the amount 
of a given type of molecule in air in the units of mass 

per unit volume, it is necessary to convert it to the 

units of molecules per unit volume to be able to use it 

in the computer program. This can be done through the 


meat Lon 


oe Na 
3 f 


M. 10> 
ale 


where oO. has the units of kilograms per cubic meter, Na 
is Avagadro's number which is 6.023 * 10 eae molecules per 
mole and M. 1s the molecular weight in grams per mole of 
the corresponding constituent. Now the amount of the oe 
constituent of air is given by the formula above in units 
of molecules per cubic meter. Final conversion may be 
made into units of molecules per square centimeter per 
kilometer which may be used directly as input to the com- 


puter program presented in Appendix A. Ws may also be 


expressed in terms of P and T as 


x. 
w. = 2.0982 x 10°? Ft 
aU Sb M. 


(3-1) 
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where P 1s alr pressure in millibars, T is temperature in 
degrees Kelvin, x and M. are as defined before. Now W, 
is the amount in molecules per square centimeter per 
kilometer for the woe type of air molecules. 

Using the equation (3-1), the amount of carbon dioxide 
which has x = 330 and M, = 44 1s 5.5943 x ia molecules 
Pee square centimeter per kilometer in the air where 
P = 1017.5 millibars and T = 286.22 degrees Kelvin. 

The densities of non uniformly mixed gases are more 
Semoltcated to determine. For ozone at sea level, it is 
taken as 6 * one grams per cubic meters or 7.5288 x one 
molecules per square centimeter per kilometer at any given 
temperature and pressure. 

The amount of water vapor and its partial pressure can 
be found from the dew point temperature or relative humidity 
in the air [{Ref. 15]. Different dew point temperatures 
and corresponding density and partial pressure of water 
vapor are tabulated in table ll. 

The water vapor density determined by interpolating 
maemetable 11 may be converted into number density (per 


Square centimeter per kilometer path) as previously 


described. 
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Although several band models were introduced by many 
workers in order to fit the molecular absorption lines in 
the atmosphere, line by line calculations are made here, 
making use of the experimentally found and compiled data. 
The molecular absorption coefficient is the most important 
meeemuUatlOn coetficient to be calculated, since mainly it 
is the one which determines the atmospheric windows which 
ere good enough for propagation of light. Also it depends 
strongly on the frequency of propagation. For any fre- 
quency of laser radiation, the molecular absorption coeffi- 
client can not be computed by interpolating between two close, 
known values. It can be computed either by the line by 
line method (infinite resolution) or by interpolation from 
calculations with resolution less than 0.01 a Therefore 
meemeacoomeed, Values GE molecular absorption coefficients 
May De used directly. If any other frequency or condition 
1S needed, the computer program presented as Appendix A may 
be used. Atmospheric windows in the 3-5 micrometer and 8-14 
micrometer regions are mainly available for atmospheric 
transmission. HF/DF laser wavelengths fall in the near 
macmmmedie iImfrared region of the spectrum. In particular 


the P.(8) transition of the DF laser with wavelength 3.8007 


5 
micrometers shows relatively high transmittance falling 


in the 3-5 micrometer region of the spectrum. The molecular 
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absorption coefficient for this particular wavelength is 


peo to4 km t 


in the midlatitude winter model where the 
molecular scattering coefficient is essentially zero. 
This shows 99 percent transmittance along a one kilometer 
path. 

Calculations intended for use in the propagation studies 
in the Monterey Bay area make use of the known range of 
conditions in Monterey, taking into account the yearly and 
monthly averages of air pressure, air temperature and dew 
point temperature in the Monterey Bay area. These yearly 
and monthly averages are taken from [Ref. 16} and [Ref. 17]; 
the quantities used in calculations are tabulated in table 
12. In this table the mean temperature is used which is 
the mean value of maximum and minimum average temperatures 
Seaeeesponaing to day/night condition. 

Also the amounts of air constituents used in the calcu- 
lations are reported in table 17. 

Table 9 shows the calculated molecular absorption for 
some selected wavelengths for the midlatitude winter model, 
using the computer program presented in Appendix A. Table 
/moncontains the molecular absorption coefficments for yearly 
averages in the Monterey Bay area. 

Tables 14, 15, 16 show the molecular scattering and 
Mpeeiemvemeat 1.06, 3.8007 and 10.591033emicrometers 
respectively for the monthly averages in the Bay area. 

Since the molecular scattering is negligible, it is assumed 


to be zero at 3.8007 and 10.591033 micrometers. 
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iitewacrosol extinction coefficients are calculated 
meee ene rural and maritime models using the size distri- 
bution given by equation (1-6a) and presented in figure 8 
and figure 9 respectively. The maritime model with 753 
Bea Spray iS normalized so that the total number of parti- 
cles is about 44 per cubic meter. These results agree with 
the previously reported values [Ref. 12] within 1.8 percent. 

The relative humidity and wind velocity dependent size 
Soe ri bution is used in figure 10 and figure 11, normalizing 
eee Gurves £O a visibility of 23 kilometers. In both 
figures 71% sea spray and 29% rural aerosols are assumed. 
Although at 0.55 micrometers, the rural aerosol refractive 
indices are outside the approximation range, the results 
seem to be valid since the percentage of rural aerosols 
ns LOW. 

As seen in figure 10, increasing wind speed increases 
the attenuation in the longer wavelength region. Similarly 
Meeraure Ll, increasing relative humidity increases the 
attenuation in the longer wavelength region. 

Figures 12 through 19 show the transmittance along a 
path length of 10 kilometers, due to the molecular absorption 
Sm ly . 

Meee ssi2, 14, 16, 18 are produced, using a spectral 
mesokmecon Or 0.01 emt for use with laser sources. 

Bedeqraded resolution is used in figures 13, 15, 17, 19 


by averaging the transmittance at each point over a spectral 
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bandwidth of 0.5 aa InemerOre theg@tagqubes=:s, [Seely 
19 are for use with the broadband sources Sil. 

mae One Of the figures 12, 13, 16, ie, 18, Po veovers 
meopeceral interval of about 25 cm™ + near the wavelength of 
Me Mrerometers (i.e., 2631.09 emss)s “The figures 14 and 
15 are for the region about the wavelength of 10.6 micrometers 
fete, 943.936 cm +). 

Figures 12 through 15 show the transmittances for 
yearly averages in the Monterey Bay area. 

Piggies 6 through 19 are produced for the different 
partial pressures of water vapor, having the same air 
temperature of 290 degrees Kelvin and the air pressure of 
moGS meltibars. The partial pressure of water vapor in 
meeutes 16 and 17 is 7.48 mbars, and it is 15.48 mbars in 
figures 18 and 19. 


Hsing cne high resolution of 0.01 cm + 


and a degraded 
resolution, a sensitivity analysis due to the air pressure, 
temperature and partial pressure of water vapor has been 
carried out near the wavelength of 3.8 micrometers. 

It has been observed that the change in the air 
pressure and/or temperature which could affect the line 
widths, line intensities and amounts of the absorbing air 
molecules, does not contribute to the transmittance signi- 
ficantly, at the region near 3.8 micrometers. The reason 


could be the small number of absorbing lines which are 


mostiy due to water vapor. Comparing figures 12, 16, 18 or 


ok 





figures 13, 17, 19 it can easily be seen that the change 
of the partial pressure of water vapor effects the trans- 
mittance Significantly. 8 millibars change in the partial 
pressure Of water vapor changes the transmittance of the 


order of 6-7 percent at 3.8 micrometers region. 
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APPENDIX C 
CeteuraATrON OF EXTINCT@ON EFFICIENCY BY USING 
THE VAN DE HULST APPROXIMATION 
This approximation was first developed by Van de Hulst 
fers 13} following the Huygens principle. It is good for 
Spherical particles of arbitrary size but [m| > l. 
Miewcsetractive index of the particle is denoted by m 


which has real and imaginary parts. 


Defining the normalized size parameter o and absorption 


parameter g such that 


18. 
Il 
I) 
{2 
a) 
; 
t- 


the extinction efficiency may be written as 


Oy = 2- 4 <o5 3 exp{-optan g] sin (0-9) 
ae 1 (SIS Sere 2g = exp [—pran =G) cos (9-2g) } 


This uncorrected extinction efficiency oe overestimates 


the cross section for small o and underestimates it alg! 
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varying degrees as 9 approaches and surpasses the value 
Pore cOrresponding to the first maximum in Ons: 
A correction to the extinction efficiency is made in 


the form 


Beene D>, the correction factor, is given by the expressions 


2 


ee) 0.2p-n_ +1 
wee a Titio) RED PS LD) 
ie!) 
_ ig eS Zz 4.08 
ee When SRS? Seiten ¢ 
eae) 
7 x 4.08 4.08 
oo én (1+3tan ee wee ep Se slia Si = © S T¥tan G 
eee (mn —1) 
eS r c 4.08 
ae n Pf (g) Belsiio2 1) ene l+tan g 
where 
2 
meee —= 1 + 4 tang + 3{tdan Gg) 


Wereman the limits of 


ra 
A 
B 

[A 
fear 
Wn 


woe 


oO 

fers 
By 

[A 
oS 
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oa meomaccurate Wlthin an error of about = 0.905 oe 


for a wide range of sizes and types of index of refraction. 
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HaSLE Lf. “MODIFICATION FACTOR FOR COz LINES 


ve= x 
L/CM 
0.0 1.00 
0.5 1.00 
0.6 0.96 
O.7 0.89 
0.8 0.32 
0.9 Clam 
1.0 0.70 
eZ 0.50 
1.5 0.50 
2-0 Q. 41 
Zoo 0.34 
aie 0 0.31 
Die 0.29 
840 0.23 
10.0 0.19 

45.0 0-00 
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TaBEE 2. 


wee eCULE 


H20 
CO2 
03 
N20 
CO 
CH4 
02 


J 


5 1.0000 
ee 1. 0095 
4.5 1.0040 
eo) be Oi.1 0 
0 1.0000 
iee> 2 «90000 
C 2.0000 


TEMPERATURE 


1.0000 
Ae Oe 
ae. 007 © 
1.0300 
1.0000 
1.0000 
1.0000 


1.20000 
1.0327 
1.0150 
1.0480 
1.0000 
1.0020 
420009 
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VIBRATIONAL PARTITION FUNCTIONS 


(DEGREES KELVIN) 


Zo 


ke 0000 
1.90502 
beo22u 
i. 07 20 
420000 
ib OO 16 
1.0000 


1.0000 
Eau? ig 
1950339 
1.1000 
1.0000 
1.0040 
1.0000 


2c 


1.000 
eo oS 
1.046 
Vode { 
1.000 
1.007 
1.000 


325 


1.0000 
1.1269 
1.0660 
Loa hee 
1.0000 
1.0110 
1.90010 





TABLE 3. SELF BROADENING ABSORPTION COEFFICIENTS FOR 
WATER VAPOR IN THE 3.5-4.2 MICROMETER REGION 


v C(v,296) 
1/ CM L/(PR -CMIATYM. 
Zo a0 0.230 
2400 0.187 
2450 0.147 
2500 O.117 
2550 0.097 
2600 0.087 
2650 0-100 
2700 ei 20 
2750 0.147 
2800 0.174 
Zo 0 0-200 
2900 0.240 
2950 0.280 
32000 G2530 
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1:6} cappegegts 





TABLE 4. CONCENTRATIONS OF UNIFORMLY MIXED GASES 
BN LHe AER 


CONST ET Gent CONCENTRATION 
POMPer SVELUME 


Gaye 2.50 0E +2 
Wes 2-300E~1 
oe (GOL =2 
CHS 1 .oU00E +0 
02 OOS) See 2. 
N2 Te ey OSE 3 bs) 


iat 





ZABLE 5. 


MODEL 


TROPICAL 


MIOLATITUDE 


SUMMER 


MIOLAT ITUDE 
WINTER 


SUBARCTIC 
SUMMER 


Sus aRcric 
WINTER 


UeSe STANDARD 


Bike ore 


MB e 


L073 


1013 


YOu 


One 


L013 


OMS 


TEMP. 


DEG. 


ele 


300.0 


294.0 


Zeer 


Cones 


283.1 


TZ 


DENSITY 
GR/M3 


57 


Phen 


i220 


Loa 2 


MODELS OF ATMOSPHERE AT SEA LEVEL 


AZ 


GR/M3 


CaO! 


14.9 


Be 
GR/M3 


5160 E=-5 


6&.Q§=5 


4.SE~5 


4 ol ote) 


5 .4E~5 





MeBLE 6. SIZE DISTRIBUTION PARAMETERS OF AEROSOL MODELS 


RURAL Oe SoCo 5 io. O0OSU 0.415 25 52—-ce Oe U 6 04415 
(CONTINENTAL) 

MARITIME 1 Oo3u Oo% 
OSee SPRAY) 


Meer oe NORMALIZED TO 2 PARTICLE/CM3 
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meine: /. Hr LASER FREQUENCIES [10] 


VE BRATIGQNAL ROTATIONAL EML Ss TON 
TRANSITION ID TRONS TEIENSIO FREQUENCY (1/C™) 


1-0 Pb 3652.50 
PT 3 ganas 
pg 3553.80 
p9 3542.20 
P10 3489.59 
P11 3436.12 
P12 3381.50 
P13 3326.21 
P14 3269.90 
P45 3212.80 
BET p2 3708.86 
p 3 Beene 
Pe AAO 7 
Ps Si aa 
Pé 3531.31 
P7 SAR GG2 
pg 2425.17 
ps aeaeeee 
P10 EE a 
Pil 3282.86 
P12 3230.18 
P13 3176.60 
P14 Biles is 
PLS 3067.22 
ve p2 ENA, Sih 
p3 3503.80 
P 4 ES yy Be 
Ps 3418.16 
Pé BTEC 
27 3327.73 
Pa 3280.4 
ae Ps DOL. 32 
Dé 3219.50 
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TABLE 7. HF LASER FREQUENCIES [10] (Continued) 


PT 3175.34 
P8 3130.09 
pg 3083.83 
5-4 P4 3150.67 
Ps Sige. 34 
P& 3068.63 
PT 3026.21 
Pg 2982.51 
Pg 2937.79 
5—5 P5 2961.68 
PE 2921.74 
P7 2380.70 
P83 2838.59 


V2 





TABLE 8. DF LASER FREQUENCIES [4] 


VIERAT IONAL ROTATIONAL EMISSION 
TRANSITION [D0 TRANSITICN ID FREQUENCY (1/CM) 


aie Pi 2884. 934 
P2 2862.652 
P3 2835.779 
P4 ZElGe. 362 
P5 2792.437 
P& Z27E7.¢14 
P7 2743.028 
PR 24s 526 
P9 26591.409 
PLO 260652200 
Ea 2638.36 
Pili ZO it~ei2> 
P10 2584.910 
P14 2557.00 
Eto 2527-060 
PL6 2498.020 

atk p32 EL TE Or Bia, 
P4 2427-330 
2H) 2703.980 
P6 2680.280 
P? Zesss oO 
PS Come © 70 
Pg 2605.870 
Pi0 25605160 
Pil 2553.970 
Pll 25276479 
P12 25002320 
Pié 2417.270 

o-e Js) Zoo2e 170 
P4 240.040 
ps 2617.410 
PE 25942230 
P7 257 Oe 0 
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TABLE 8. DF LASER FREQUENCIES [4] (Continued) 


P8 2546.370 
po 2521.810 
P10 2496.610 
Pil 2471 .340 
P12 | 2445.290 
P13 2415.020 
P14 2392.460 
4-3 PS 2532.500 
Pé 25094860 
P7 2486. 830 
P8 24634250 
pe 2439<«290 
Pic 2414-890 
5-4 P7 2404.630 
7-6 Pg 2222.680 
Pi0 2177.990 
PLL 2155.030 
P12 2121.680 
8-7 P7 2165.920 
PB 21442300 
Pe 2123.240 
P10 Zi .2 70 
P1I 20264140 
Pi0 2032.010 
5-8 P6 2108.480 
P7 2088.340 
P8 20672760 
P10 20254360 
Pil 2003.560 
P12 1981.330 


dl 





TABLE 9. MIDLATITUDE WINTER MODEL MOLECULAR ABSORPTION 
CORP ETCTENTS FOR INFINITE RESOLUTION 


FREQUENCY MOLECULAR ABSORPTION 
(1/CM) COEFFICIENT (1/K™) 
924.975 4276-1 
942 .384 6.17E-2 
944.195 6.61E-2 
945.981 6.56E-2 
947 .743 7.22E-2 
959.393 3.456-2 
973.289 7649E-2 
974.623 6 .96E-2 
Besos! 1.32E-1 
SCTE Pals 7.28&-2 
1046.854 8.27E&-2 
1048.661 8.97E-2 
1050.441 $.81E =2 
1075.988 8.902-2 
1077.303 8.21E-2 
Hove 6591 % 375-2 
2177.990 2.389E-2 
2414.890 1. LGE-2 
2463.250 1 .26E-2 
261.125 7.19E-2 
2617.410 7.1 9E-3 
2631.090 1.04E-2 
2638.396 7.30E-2 
2640 .040 1.67E-2 
2 655.970 2.24E-2 
2662.170 1.246-2 
2665.200 1 -64E-¢ 
2703.980 9.55E-3 
Aah B ake 2.37E-2 
27 27.380 1.386-2 
2743.028 1166-2 
2750 .050 yeoLe—2 
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meer 10. 


WAVELENGTH 


(um) 


025145 
0.55 

Oees23 
0. E543 


Riv hAeenve INDEEGES OF AEROSOL PARTICLES 


NST nes SOLUBLE DUST LEEKE 


Gar? O 
0.030 
0.008 
0.008 
0.008 
0.003 
0.003 
0.008 
0.008 
0.008 
0.008 
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0.098 
0.008 
0-908 
0.008 
0.009 
Ye O09 
0 .Ui3 
Org 
0.040 
0.013 
Cot 
0.011 
0.012 
0.014% 
0.016 
0.92% 
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0.039 
0.042 
0.055 


See] veIKe 


42 5Q 
ie 2 
1.74 
et 
(i 0S 
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Veto 
eto 
PA c=) 
eS 
(oo HS 
oie 
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4ef9 
1.8090 
Dayo 
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ecb of 
6 a e 
Ve) 
on Iv On £ 


he 
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Oy —O® © ww so 
ol 


ht fo 
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Ona 


C.45 
O+? 
C.47 


Gare. 


0.+5 
0 45 
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Oe+3 
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0.44 


C.45 | 


Oe+6 
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Jose 
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PapbGE 10. 


169 


Sez 


REPRAGTIVE INDICES CF AEROSOL PARTICLES 


[Ref. 
L.200 


1. 010 
1.300 


1 


C.065 
0.100 
Gece 
0.290 
O. 370 
06420 
Geico 
Gso> > 
0.090 
G.070 
0.050 
0.047 
0.055 
G02 2 
0.073 
C.i90 
0.200 
0.4150 
0.240 
0.180 
0.170 
0.220 
0.230 
0.240 
Ce 280 
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G.500 
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O. 590 


0. O+0 
0.074 
0.090 

ene 
Gel40 
0.150 
0.162 
0.162 
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02420 
Oxrvi> 
0.100 
0.090 
0.iL00 
0.085 
0.100 
0-100 
0.100 
0.100 
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O20 
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Oe 240 
0.2329 
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0-600 2 
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fap tl. WATER VAPOR PRESSURE AND DENSITY FOR DIFFERENT 
DEW POINT TEMPERATURES 


DEW POINT TEMP. VAPGR PRESSURE VAPOR DENSITY 


(DEG. KEL.) (MILI BARS) (GRAMS/M3)} 
229 0.189 0.176 
238 0.314 0.286 
243 o- 00 On 452 
248 0.807 0-705 
253 1.254 1.074 
258 Peo igo O05 
263 22863 Paty ays 
225 4.215 3.407 
fe 6-108 4.347 
278 Sie tieg 0 Tet! 
23 3 LO Saeco 
288 172044 Mees o0 
Ue) © 2re 5) 1 17.500 
253 Sere ts 23-050 
B03 426430 5302380 
308 56 236 Spe) ie) G0, 
31> poe rete Seek 20 
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TABLE 13. MOLECULAR ABSORPTION AND SCATTERING COEFFICIENTS 


FOR YEARLY AVERAGES IN THE MONTEREY BAY AREA 


FREQUENCY TRANSITION MOL. ABS. MOL. SCAT. TRANS. ALONG 
ID C2 eit. @UEr. 1 KM. PATH 
(1 /CM) (1/KM) (i /KM) (3) 
944.155 P(20) C02 0.1520 = 25.36 
goees eS P{22) COZ 0.1460 = 86.44 
2471.2340 P{il) OF 0.0230 = G1w 72 
24966510 P(1L0) ODF 0.03 84¢ = $8.17 
Pez iawtiO P(i2) OF 0.0172 = 2G a) 
Bose stO P(t) = 6 6(DF 0. 0129 = S/S) SOE 
BesOetoO P(10) OF 0.0293 = Cet 
Boua6eotT) P59) DF 0.0240 = SoG 
woowe09O P( 8) DF O. 0228 = BU tele 
Meoce 396 P{(IT1L) OF Oe Zee = GO el / 
Moose 97 CG PCT) DF OE) Syl = 94.33 
moooesoc P(10) £4xOF CAS) aie = $6.49 
moouec SC PP (6) DF 0. G547 = 4.63 
ese 260 8 P(5) Big 0.0208 = $729% 
Meee sS Pd) DF 0.0930 = Ere te 
mee I2S PCs) DF 0.0246 = ies og | 
2f767.914 Pd) DF 0.0594 = $4 223 
ero eetS i =6P( 5) DF 0.0455 = Ser Tile 
3373-460 P(5) HF 0.1240 = 88.37 
Seo t oP (3) HF Ore = 70.72 
3483 .63C PT) HF 25 30551610) = 1.32 
Beeewore) PLS) rae 8625000 = 0-00 
Berheti OC P(5} ‘at 35.0000 = 9.00 
BoooeoO0 P(8) HF 67.2000 = 0.00 
3622-710 P(4) HF 3526000 3 0.00 
9433.3962 0.000014 0.0008263 ©9.$2 





TaBboe 14. 


JANUARY 
eon AR Y 
MARCH 
APRIL 
MAY 

JUNE 
JULY 
AUGUST 
Ser , EMBER 
SE TOSER 
NOVEMBER 
BEeCeMBeER 
=a RLY 


MOLECULAR ABSORPTION AND SCATTERING COEFFICIENT 
Pie Oo MICROMETERS IN THE MONTEREY BAY AREA 


MOL. ABS. 


COEF. 


(1/KM) 


ea Claws 
9 ./0Q0E-6 
elit Case 
Pe 7 OOE=s 
iil gets (8 (CNS— 7s 
13.300E-6 
Mom 700 Es 
eG. 5006-5 
14-100E-6 
iz sO OE=6 
19.400 6 
9 41LOE-6 
12.4Q0&-6 


MOL. SCAT. 


Coe 
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Se S49E=4¢ 
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8e27ZE—-4 
So 223E=4 
Be cs2t= 
S.220e-4 
Seki cease 
Bee 50e=4 
&@.304E—4 
8.344E-4 
e272 8SE=4 
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TOTAL MOL. 


EXTINCTION 


(1 /KM) 


8. 48EE-% 
8-474E-4 
8-4465-4 
8-42 7E-4 
B.2 30-4 
Se361E-4 
8e3 7LE=4 
S19 6 DE +4 
8.-313E-4 
Sel oL eas 
8 e408E-4 
8. 428E-4 
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Se ae 
Bele 
S20 Se 
oe 2 
ee 
Bela Ss 2 
S592 
mea? 
ope) (Sli 
She a) 
SS OS 
Boe 2 
99.92 
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Peewee to. MOLECULAR ABSORPTION COEFFICIENTS AT 3.8007 


MIGROME TERS “EN THE MONTEREY BAY AREA 


MOL. ABSe TRANSe ALONG ORDER OF 
COEF. ewe eo cast MON 
(1/KM) (2) TG PROPAGATE 
JANUARY 0.0194 93.08 2 
FEBRUARY 0.0199 936 02 “ 
MARCH 0.0199 7oa05 3 
APRIL OS Ome Si .38 6 
MAY OsOce2 Siaig 7 
JUNE 0.0242 OF ok e 
JULY 0.0249 I7¢ 54 ee 
AUGUST 0.0256 97 6438 a2 
SSP 1) Baal See 0.0247 Phe sre 0, 
OCTOSER 0.022 7 Shots 8 
NOVEMBER 0.0202 33.00 2 
OeGe sex 0.0153 98.09 i 
YEARLY 93.23% 


0.0218 
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feito. MOLECULAR ABSORPTION COEFFICIENTS AT 10.591033 
MICROMETERS IN THE MONTEREY BAY AREA 
MOL. ABS. TRANS. ALONG ORDER OF 

CUERT 2 KM. PATH S8EST MONTH 
(2/KM) (3) TI PROPAGATE 

JANUARY Oei2e2 33.52 Z 

Peon UA Oe age: 83. OL 4 

MARCH 0.128 S36 U2 3 

APRIL 0.147 Sold 6 

MAY 0.359 32 aes) if 

JUNE 00.135 836405 2 

JULY 0.195 S2.52 Ha 

AUGUST 0.204 31.54 Aaa 

SEPT eVGeR 0.193 8243 10 

GiGy Geer O. 26> 84.78 8 

NOVEMBER Oars i eero = > 

Mecers ER Cite 83.57 AL 

YEARLY Oe ers 35.36 
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Teapots. AmROSOL EXTINCTION COEFFICIENTS FOR MARITIME 
MODEL WITH 75% SEA SPRAY AND 25% RURAL AEROSOLS 


WAVELENGTH AEROSOL EXTINCTION AEROSOL EXTINCTION 
(MICROMETERS ) COEFFICIENT (1/KM) COEFFICIENT (1/KM) 
CALC UMA TED REPORTED [Ref. 12] 
2A, 0.11330 Ooi Seve 
255 0.09895 0. 09962 
Ere 0. 08258 
3.0 0.10430 0.10426 
3.2 0.10800 
3.3923 Oe G20 
ols 0.09945 0.09899 
3.75 0.09269 0.09193 
4.0 0.08802 0.08670 
4.5 0.07875 
565 0.05007 0.05928 
6.0 Genes 5515 0.05485 
6.5 0.05825 
fee 0.04304 0.04758 
as 0.04205 0. 04063 
Sar 0.03958 0.03960 
a.5 0.03572 0. 04045 
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Poca. MODEL OF AN ELEMENTAL SCATTERER 
b. CREATION OF AN INDUCED DIPOLE MOMENT 
BYigaN biG TRG Ger tei 
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